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ABSTRACT 

Two aspects of filamentary molecular cloud evolution are addressed: (1) Exploring analytically 
the role of the environment for the evolution of filaments demonstrates that considering them in 
isolation (i.e. just addressing the fragmentation stability) will result in unphysical conclusions about 
the filament's properties. Accretion can also explain the observed decorrelation between FWHM and 
peak column density. (2) Free-fall accretion onto finite filaments can lead to the characteristic "fans" 
of infrared-dark clouds around star- forming regions. The fans may form due to tidal forces mostly 
arising at the ends of the filaments, consistent with numerical models and earlier analytical studies. 

Subject headings: methods: analytical — stars: formation — ISM: clouds — gravitation — MHD 



1. MOTIVATION 

Recent observations have re-emphasized the promi- 
nence of "filamentary" (at minimu m, non-spherical) 
structures in molecu lar clouds ([A ndre et al] 120101 : 
Moli nari et~aTl 120101 : lArzoumanian et al.l l2011f ). 
While certainly a con f irma tion of earlier findings 
([Schneider fc Elmegreenl Il979h . higher observational 
resolution probes successively deeper into the structure, 
and thus into the physical parameter space occupied by 
these "filaments". 

Filaments are the most beneficial geometry for lo- 
cal gravitational collapse to occur before global collapse 
(|Pon et al.l l2011f ). The global free-fall time increases 
with respect to t he (spherical) three-dimensional case 
(jToala et al.l 120121 ). and - equivalently - the ratio of the 
local over the glob al collapse time scale decreases with 
reduced dimension (|Pon et al.ll2011f ). 

Self-gravitating molecular cloud s rarely exist in com- 
plete isolation. iHeitsch et al.l (|2009af ) showed for 
model clouds of s imilar parameters as the Pipe Nebula 
(|Lada et al.l 12008) that once the clouds are not consid- 
ered in isolation, the gravitational potential generated by 
the cloud mass will lead to accretion. In fact, Pipe-like 
clouds are likely to accrete their own mass with i n a dy- 
namical (read "free-fall") time. IHeitsch et all (j2009af ) 
also demonstrate that while their models show clear 
infall and accretion, the kinematic signatures of such 
motions can be easily hidden if the medium is struc- 
tured, or by projection effects. Observed morphologies 
and ki nematics su ggest gas infall, from the filament-core 
level ([Kirk et al.l 120131 see a lso e.g. Hacar & Tafalla 



2011 ), on to low-mass clouds (Ballest eros-Paredes et al. 



1999; iPalmeirim et al.l 120131 for Taurus) up to clouds 



harbouring massive star- form a tion regions dMotte et al.l 
20051: [Sc hneider et al.l 120101 : iHennemann et al.l l2Q12t 



Pitann et al.ll2013f ). These studies also suggest that the 



prominent filaments are hosting most of the star forma- 
tion, and thus geometric properties of filaments could 
cont rol the star fo r matio n process to a large extent. 
lArzou manian etall ()2011f ) show that for their sample of 
filaments, the filament width does not depend on the cen- 
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tral column density, and speculate that this could be an 
effect due to accretion. Gas accretion along filaments has 
also b een predicted based on simulations (Bal sara et al.l 
l200lh . 

The notion that molecular clouds do not exist in iso- 
lation, and thus will feed on their environment, raises 
two issues that we explore for a better understanding of 
molecular cloud evolution. 

(1) The Fate of an Accreting Filament — Gas accretion 
onto a filament will not only grow the cloud's mass and 
thus reduce the fragmentatio n timescale, but it a l so wi ll 
drive internal turbulence (Klesse n fc Hennebel le 2010). 
The relevance of accretion over fragmentation will be ex- 
plored in $21 Already a comparison of the characteris- 
tic timescales shows that accretion ought not to be ne- 
glected. 

(2) Free-fall Accretion onto a Filament — Accretion onto 
a filament will lead to characteristic kinematic signa- 
tures that will be explored in $3] Specifically, tidal 
forces around a massive filament can cause "gravitational 
streamers" , i.e. fan-like structures converging at the ends 
of a filament. 

The considerations presented here are purely analyti- 
cal, and they are somewhat speculative. The goal is to 
explore an extreme view, namely that free-fall accretion 
governs the evolution of molecular clouds. $4] highlights 
some of the arguments for and against this assumption. 

2. EVOLUTION OF AN ACCRETING FILAMENT 

In this model, the evolution of an accreting filament 
is determined by two timescales: the gravitational frag- 
mentation timescale along the axis of the filament, and 
the accretion timescale, i.e. the timescale it takes the 
filament to reach line masses 
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above the critical line m ass m cr (e.g. [ Ostriker 1964]), or 
- in the nomenclature of iFischera fc M artin (2012) - to 
reach 
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We will denote the isothermal sound speed by c s 
the line mass by m. 



and 



2.1. 



Comparison of Instantaneous Timescales 

The most unstable scale of an isothermal, hydrostatic 
cylinder of infinite length is given by 
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with the central density p c = [iVfiHin c . The cylinder will 
fragment at this scale on a timescale of 
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(iNaga sawa 1987; Tomisaka 1995), d ue to a varicose or 
sausage instability (see Uackson et al.l[2QlQ, for an appli- 
cation to the "Nessie" filament; also M4.1.2p . 

The free-fall accretion timescale onto a filament can be 
derived by the ratio of the line mass m over the accretion 
rate rh through a cylinder circumference 2ttR, 
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Following the discussion of He itsch et al.l ()2QQ9ai see also 
Palmeirim et al. 2013) to determine the accretion veloc- 
ity vr, the gravitational acceleration of a gas parcel in 
free-fall onto a cylinder with line mass m is given by 
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Assuming a steady-state infall, this yields a radial veloc- 
ity profile via 
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Four parameters remain to be determined: the parcel 
density p, the line mass m, and the radii R re f and R. The 
ambient gas density p ext will serve as the parcel density, 
assuming an average over a highly structured medium. 
The line mass m is set to the critical value for now - 
awaiting a more detailed treatment in the next section. 
The choice of the radii is trickier. For now, suffice it to 
say that vr depends only weakly on R re f/R, so that this 
choice is of lesser importance. The integration constant 
R re f can be identified with the starting position of the 
fluid parcel, located at - somewhat arbitrarily - R re f = 
2 pc. The radius R itself is set to a characteristic filament 
widt h from observations, e.g. a few tenths of a parsec 
(e.g. lArzoumanian et aLll2Qll[ ). 



Combining all expressions leads to 
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Thus, the accretion timescale (eq. [9j) and the fragmen- 
tation timescale (eq. H]) are of the same order. Fila- 
ments that do not exist in isolation accrete gas on simi- 
lar timescales under which they fragment. Yet, the above 
estimate suffers from the somewhat random choice of pa- 
rameters, begging a more detailed discussion. 



2.2. Time- dependent Timescales 

The next step beyond considering instantaneous 
timescales involves integrating the expression for the 
mass accretion rate using equation [8j 
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Thus, the accretion rate will depend on the evolutionary 
stage of the filament via the line mass m(t), and on the 
surface radius R. The latter is set to the filament radius. 
This would be formally infinite for an isothermal, hydro- 
static cylinder. Yet since filaments do not exist in vacuo, 
it seems reasonable to define the (outer) filament radius 
Rf as the radius R at which the radial density profile 



p(R) 



has dropped to the ambient density p ext , i.e. 
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For an isothermal cylinder, p = 4 . The charact eris- 
tic ra dius Rq (equivalent to Rfi a t in lArzoumanian et al.l 



(120111) ) is set to the isothermal core radius (jOstrike: 
196|), 
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using equation [TJ Note that t his expre s sion d iffers by a 
factor of 8 from that given by 10 striker! (|1964f ) . 

Unlike in the analysi s of pressurized, hydros tatic 
isothermal cylinders by Fisc hera fc Martini (J2012), we 
specify the steepness of the filament profile explicitly by 
setting p, instead of relying on an external pressure to 
flatten the profile. We discuss the rationale and physical 
interpration in §4.1.11 

The line mass m(t) can be found by integrating over 
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the density profile and by using equations [12] and [13j 
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Note that m(Rf) > for all values of p > 0, but that m 
will not converge for p < 2. Yet, this is of minor concern 
here, since the radius Rf will always be finite by con- 
struction (see eq. [12]). For p = 4, m(Rf —> oo) — » m cr . 
The initial condition for equation [10] is set by assuming 
an initial central density p c , and a (constant) ambient 
density. The central density as a function of m is given 
by 
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for p = 2, 
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where we used equation [2] 

Figure [1] (black lines) summarizes the results of the in- 
tegration for p = 3 (i.e. flatter tha n isothermal, and 
steep er than values o bserved by [A rzouman Tan et al.l 
2011, see, however, lHacar fc T afalla 2 Oil) regarding 
evidence for isothermal profiles), R re f — 2pc, no = 
100 cm -3 , an initial central density of n c (0) = 200 cm -3 , 
and for T = 10 K. For this parameter choice, the filament 
reaches criticality after « 1 Myr (Fig. [TJi). Accretion 
dominates the filament evolution until shortly after that 
time (Fig. [lb, black lines). Other choices of p give similar 
results. 

The core radius Rq (Fig. [l]i, dashed lines) decreases 
monotonically, because of the growing central density 
(Fig. [!£). The filament radius Rf (defined as the radius 
at which the density profile drops to the background den- 
sity) increases to a maximum, and then drops, since for 

large p c , Rf ex Pc 6 , and for small p c > p ext , Rf ex pi . 
Also shown in Figure QJ1 is the most unstable wave 
mode A max , normalized to its initial value of « 4 pc. 
With increasing central density, \ m ax decreases. 

2.3. Effects of Accretion- driven Turbulence 

iKlessen fc Hennebelld (|2010l ) argue that in many as- 
trophysical objects turbulence can be driven by accre- 
tion. For molecular clouds, this point has been made by 
various gro ups based on simulations of flow-driven cloud 
form ation ([Vazquez- Semadeni et al.ll2007l : iHeitsch et al.l 
2008, see also iField et all (j2QQ8h for a more systematic 
approach). In this scenario, turbulence is a consequence 
of the formation process initially, while at later stages, 
global gravitational accelerations drive "turbulent" mo- 
tions. 

IKlessen fc Hennebelld ((2010) estimate the level of tur- 
bulence driven by accretion (their eqs. 2, 3, and 23). For 
the purposes here, the characteristic length scale is 2Rf, 
the accretion velocity v(Rf), and the driving efficiency 
e = 0.1 (see their eq. (23)). Then, the turbulent velocity 
dispersion is given by 
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Fig. 1. — Filament properties against time for the isothermal 
(black) and the turbulent (red) case, (a) Line mass over critical 
value (eq. [TJ. (b) Fragmentation timescale Tf (solid line) and ac- 
cretion timescale r a (dashed lin e) . (c) Central filament density n c . 
(d) The filament radius Rf (eq. I12|) and the core radius Ro (eq. fT3|) 
in parsec, and the wavelength of the most unstable mode (eq. [3j) 
normalized to its initial value of 4.5 pc. (e) accretion velocity v(Rf) 
(dashed line) and t urbu lent velocity a driven by accretion (solid 
line; see eq. 1161 and fl2. 3|) . 

The choice of e = .1 er rs on the generous side - 
IKlessen fc Hennebelld (|2010l ) quote values of a few per- 
cent. 

Figure []J (black lines) shows the velocity dispersion a 
against time (solid line), and the accretion velocity at 
Rf for the filament discussed above (dashed line). The 
velocity dispersion a is a factor of a few larger than the 
sound speed at T = 10K used to evaluate Ro. For e = 
0.05, a is on the order of the sound speed. 
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Yet, the above estimate of internal motions driven by 
accretion does not account for the effect of an increased 
a on the critical line mass m cr and on the core radius Rq. 
The latter will affect the line mass m of the filament, and 
thus the accretion rate. Since the internal motions are 
noticeably driven by infall at a factor of a few, a closer 
look at this backreaction seems in place. That the above 
discussion is inconsistent because of the missing effect 
of g on the filament growth is also demonstrated by the 
shape of a(t) (solid black line in Fig.QJ): just from virial 
considerations, g should increase with m. 

Replacing the sound speed used to calculate Rn by the 
velocity dispersion g has several consequence^- Since 
Ro is needed to determine not only the critical line mass 
but also the line mass of the accreting filament (eq. fl4|) , 
the accretion rate (eq. [10]) will now depend on a, and 
thus on itself. A root finder can be used to determine 
ra(cr), and thus the accretion rate m. 

As the red lines in Figure QJi show, the ratio m/m cr 
increases now at approximately half the rate, while the 
central density grows even more slowly. After 4 Myr, 
the filament has reached n c = 10 4 cm -3 , as compared 
to ~ 1.8 Myr for the isothermal case. Also note that 
the isothermal n c begins to diverge, while the turbulent 
n c does not. This is a direct consequence of continued 
accretion, driving internal motions and fragmenting the 
filament, as can be seen from Figure QJ: g keeps rising 
for the turbulent case (red lines), driven by continued 
infall due to the ever-deepening gravitational potential. 

The fragmentation timescale r/ (Fig. [TJd) and the most 
unstable wavelength (Fig. [l]i) depend on g only through 
the central filament density n c . For the fragmentation 
timescale, this is obvious from equation 2J and for the 
most unstable mode (eq.[3j), we use the isothermal sound 
speed, not the velocity dispersion a, consistent with our 
interpretation of turbulence not being able to provide 
pressure support. The filament and the core radius Rf 
and Ro are roughly a factor of 2 larger compared to the 
isothermal case. The filament radius does not change 
significantly, while the core radius drops by a factor of 
~5 to a few tenths of a parsec. 

Summarizing, accretion-driven turbulence reduces the 
growth rates in all quantities, but it does not eventually 
stabilize the filament: the fragmentation timescale still 
wins over the accretion timescale. Replacing the sound 
speed by the velocity dispersion g in X m ax would increase 
the fragmentation scale at a given time, but this seems 
an improper thing to do. 

2.4. Accretion of Magnetic Fields 

Fie ge fc Pudritzl (j2000al lbl) discussed in great detail 
equilibrium configurations of cylinders with a variety of 
magnetic field geometries. Here, we are interested less in 

2 This is often referred to as "turbulent support" . As has been 
pointed out repeatedly, supersonic isothermal turbulence cannot 
support a gravitationally unstable region already for the s imple 
reaso n that it leads to local fragmentation first (Klesse n et al.l 
[20001 ; [Heitsch et al.| [2001, although [Clark fe B onnell (2005) point 
out that the earlier picture of a local reduction of the Jeans mass 
may be incorrect). A perfect example of this effe ct can be seen i n 
the simulations of turbulent, isolated clouds by Clark et al. (2005). 
What is commonly referred to as "support" is in fact the effect of 
redistributing the densities in a region, leading to low volume filling 
factors and - eventually - thus to lowered star formation efficien- 
cies. 



stable configurations, but to what extent magnetic fields 
affect the growth of the filament. 

The effect of magnetic fields can be included approxi- 
mately by modifying the sound speed 
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Because of flux- freezing, a scaling of the magnetic field 
strength with density is assumed, 
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which will depend on the field geometry through the ex- 
ponent s. For fields along the axis of the filament and 
for toroidal fields, s = 1, with 
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For a uniform field perpendicular to the filament axis, 
5 = 1/2 under mass and flux conservation, resulting in 
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In this case, the magnetic pressure will stay constant, at 
its initial level. Both cases s = 1/2,1 will be explored. 
As in the turbulent case, eq.[T71can be solved numerically, 
since the RHS depends on c s via the central density n c . 

The weak scaling (s = 1/2, Fig. [2]) leads to qualita- 
tively similar results as the isothermal case - the mag- 
netic curves have the same shape as the isothermal ones. 
The only effect of the magnetic field in this case is to 
slow down the accretion, similar to the turbulent case. 
Fields parallel (s = 1, Fig. [3]) to the filament, however, 
affect the accretion properties qualitatively. Both weak 
and strong fields (/?o = 10,1,0.3) can slow down accre- 
tion, indicated by the flattening density and line mass 
curves at late times. The fragmentation timescale de- 
pends on the field strength through the central density 
n c . The length scale of maximum growth (eq.[3j) depends 
directly on the modified sound speed (eq. [T7j) - consis- 
tent with the notion that sufficiently strong magnetic 
fields can suppress gravitational fragmentation^]. 

While magnetic fields can affect the mass accretion 
onto the filament, they do not prevent it (at least 
not for reasonable choices of magnetization). Frag- 
mentation wins over accretion once the filament be- 
comes critical. We forego the discussion of turbulence 
in combination with magnetic fields. There is numer- 
ical and analytical evidence that turbulence combined 
with magnetic flux loss mechanisms (ambipolar drift 
and reconnection) efficiently red uce the d ynamic al im- 
portance of magnetic fields JL azarian & Vishniac 1999; 
Santos-Lima et al. 2010; Kim & Diamond 2002; Zweibel 



2002; Hei tsch et all 12004 



3 A different question is whether such strong fields can be reached 
in gravitationally contracting regions. 
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Fig. 2. — Filament properties against time for the isothermal case 
(black), and the magnetic cases at (3o = 10 (red), /3o = 1 (green), 
and (3o = 0.3 (blue), for a field perpendicular to the filament (s = 
1/2). (a) Line mass over critical value (eq. [1}. (b) Fragmentation 
timescale tj (solid line) and accretion timescale r a (dashed line), 
(c) Central filament density n c . (d) The filament radius Rf feq. U2|) 
and the core radius Ro (eq. [T3]) in parsec, and the wavelength of 
the most unstable mode (eq. [3} normalized to its initial value, (e) 
accretion velocity v(Rf) (dashed line) and turbulent velocity a 
driven by accretion (solid line; see eq. 1161 and fl2. 3|) . 



EFFECT OF TIDAL FORCES ON FILAMENT 
ACCRETION 



iMversI (J2QQ9L 12011( 1 discussed a detailed model of 
star formation regions of the hub- filament morphol- 
ogy (see a list of examples in IMversI I2QQ9L §2). He 
then explored a model in which clumpy molecular cloud 
gas is compressed into a vertically self-gravitating and 
density- modulated layer, preserving filamentary struc- 
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Fig. 3. — Filament properties against time for the isothermal case 
(black), and the magnetic cases at /?o = 10 (red), /5o = 1 (green), 
and /3o = 0.3 (blue), for a field parallel to the filament (s = 1). (a) 
Line mass over critical value (eq. U>. (b) Fragmentation timescale 
Tf (solid line) and accretion timescale r a (dashed line), (c) Central 
filament density n c . (d) The filament radius Rf feci. I12j) and the 
core radius Rq (eq. fT3|) in parsec, and the wavelength of the most 
unstable mode (eq. [3j) normalized to its initial value, (e) accretion 
velocity v(Rf) (dashed line) and turbulent velocity a driven by 
accretion (solid line; see eq. 1161 and £|2.3j) . 



ture in equilibrium and in collapse. iSchneider et al.l 
(2010) observed "hub-type" filaments around the mas- 
sive star- forming filament DR21, including i nfall /collapse 
signat ures onto the massive cores. ISchneider et al.l 
(2012) concluded from their analysis of the Rosette fil- 
aments that clusters form at the position of filament 
mergers or crossings, c onsistent with numerica l sim- 
ulations (iLi et al.l 12004 IVazquez- Semade ni et al.l 120071 : 
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Fig. 4. — Column density projection of 3D isothermal collapse of 
elliptical molecular cloud. Arrows indicate velocity vectors. 

ICarroll-Nellenback et al.ll2Q13h . 

As an a lternat i ve int e rpreta tion to the equilibrium 
models of iMversl (J2QQ9L 12011( 1 on the one hand and 
as an extreme interp retation of the turbulent view 
([Schneider et al.l 12012) on the other we assume that the 
gas around the massive filament is in free-fall, since the 
massive star formation region at the center must domi- 
nate the gravitational energy budget of the system. The 
" fanning-out" filaments (or "gravitational streamers") 
may be formed by the stretching of over-dense gas parcels 
( "clumpy" cloud structure) due to the tidal forces of the 
central filament. This interpretation is largely motivated 
by the notion of "focal" p oints in finite, non-spherica l, 
self-gravitating structures ([Burkert fc Hartmannl 12004). 
i.e. by the fact that gravitational accelerations peak at 
the ends of filaments or ellipses. These locally strongly 
varying accelerations not only lead to loca l com pres- 
sions due to global gravity (Bastien] ll983l ; |Li 2001), but 
they also affect the ratio of the local to global collapse 
time, rendering structures of lower dimensionality more 
pron e to suffer local collapse before they collapse glob- 
ally (|Pon et al.l 1201 1[ ). As a corollary, the y increase the 
global free-fall ti me of the structure (jToala et al.l 12012; 
iPonet al.|[2Q12bh . 

Figure 2] shows a column density projection (approxi- 
mately 2 pc across) of the three-dimensional isothermal 
collapse of an elliptical cloud. At this stage, the cloud 
has collapsed to a filament of M — 1O 4 M0, but ambient 
gas is still being accreted. Specifically, at the end points 
of the filament (just one side is shown in Fig. SJ), mate- 
rial is being accelerated parallel and perpendicular to the 
filament axis, leading to the "velocity fan" indicated by 
the vectors. The initial density structure of the ambient 
gas is slighly clumpy, not filamentary, yet, the resulting 
structures fanning out from the filament's end point are 
reminiscent of the hub filaments seen in observations. 

3.1. The Role of Tidal Forces 




x [pc] 

Fig. 5. — Trajectories for a fluid parcel under free-fall onto a 
filament, for n c = 10 5 cm -3 . The filament - modeled as a cylinder, 
see eq. [23]- is indicated by contour lines between —2 and Ope (only 
one quarter of the domain is shown). Black line segments show 
the orientation of the fluid parcel's long axis. Numbers give times 
along the trajectory in Myr. Each of the fluid parcels starts at rest 
(see also Fig. [6j)). Line segments parallel with the trajectories are 
hard to discern. 



The tidal accelerations on a gas parcel in free-fall to- 
wards a mass M at a given distance R are set by 
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with the size of the parcel Ar, and its distance to the 
mass, r. For a distance r = 1 pc, a size of Ar = 0.1 pc, 
and for a mass of M — 1O 4 M , this results in a tidal 
acceleration of 0.87km s _1 per 10 5 years. 

The effect of tidal forces on spherical fluid parcels can 
be estimated via a ballistic integration. The parcels are 
initially at rest, accelerating towards the filament. The 
filament is modeled as an isothermal cylinder of total 
mass M, over a given length scale. The density profile is 
set by 

n(x, y, z) = n + — ^ (23) 
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x 2 + y 2 + z 2 for \x\ > Xf 
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z 2 for \x\ < Xf. 



(24) 



Here, Xf = 2pc is the filament's half length, i.e. the 
filament extends from —Xf to Xf. The background den- 
sity no = 10 2 cm -3 , and central filament densities of 
n c = 10 4 , 10 5 cm -3 are explored. The core radius is set 
to Ro = 0.1 pc, and p = 3. Only one quadrant of the fil- 
ament is shown in Figure [5j but the forces are calculated 
using the whole filament. The fluid parcels themselves 
are pressure- less, thus they can suffer infinite compres- 
sion. 

We note the following from Figure First, all parcels 
are being stretched along their trajectory, as indicated 
by the line segments giving the orientation of the major 
axis of a parcel. Second, the parcels hit the filament 
inwards of its endpoint (as a reminder, the center of mass 
resides at x = 0, y = 0). Third, with the stretching 
comes an overall acceleration of the parcel. And finally, 
all trajectories are curved towards the filament. 



Gravitational Infall onto Filaments 



3.2. Predictions 

3.2.1. Curvature of Trajectories 

The most obvious difference between Figure [5] on the 
one hand, and observations of gravitational streamers (or 
hub-filaments) as well as gas-dynamical models on the 
other (Fig. H]) is that the ballistic trajectories curve to- 
wards the filament, while the observed and gas-dynamical 
ones tend to merge with the filament. A closer look 
at Figure [5] provides the following information: Sam- 
pling the (x, ?/)-plane with more trajectories shows that 
the trajectories in the vicinity of the black line actually 
cross. Though the fluid parcels are treated as pressure- 
less here, pressure forces would resist this compression 
perpendicular to the trajectory, preventing the crossing. 
Thus, a parcel's motion towards y = would be deceler- 
ated, leading to a reduced velocity component in y, and 
hence to an effective curving of the trajectory such that 
it merges with the filament. 

3.2.2. Gas Density along Trajectory 

Observed hub filame nts are more p ronounced closer the 
central structure (see iMversI I2Q09L for a list of hub ob- 
jects). Hence, the (column) density within the filaments 
is increasing. Thus, one would expect from the model 
that the fluid parcels are not only be stretched, but also 
compressed. The inverse of the volume compression ratio 
V(Q)/V(t), which is equivalent to the density compres- 
sion, for all trajectories is shown in Figure [6^. Compres- 
sion is highest for trajectories with a strong component 
along the filament axis, while it is less pronounced for 
trajectories perpendicular to the filament. Trajectories 
along the filament axis (i.e trajectories that terminate at 
one of the filament's ends) effectively see accelerations 
oc R~ 2 , while trajectories terminating further in are sub- 
ject to accelerations ex i? _1 , with correspondingly weaker 
tidal forces. 

3.2.3. Line- of- Sight Velocity Signatures 

Figure [6b shows the parcel velocities along all trajec- 
tories. These are consistent with the estimate in eq. [23j 
Yet, an observer would not see the trajectories, but radial 
velocities along lines of sight at a given time. Figure 
provides an "observer's view" of the velocity structure 
around the filament. To derive the spectra, lines-of-sight 
(as shown in Fig. [8j) are calculated through the filament 
and the ambient, inf ailing gas. Arrows in Fig. [8] indicate 
the viewing direction. Spectra are determined assuming 
optically thin tracers corresponding to two density ranges 
(10 2 < n < 10 3 and 10 3 < n < 10 4 ), loosely identified as 
12 CO and C 18 0. An initial parcel density n = 10 2 cm -3 
has been assumed. The sound speed in the accreting gas 
is set to c s = 0.2 km s _1 , and to 1 km s _1 in the central 
filament, mimicking the effect of accretion-driven turbu- 
lence. The emissivity is assumed to be proportional to 
the local density enhancement due to fluid parcel com- 
pression (see §3.2.2p . Emission from the infalling gas 
and from the central filament contribute to the spectra, 
providing an estimate for the visibility of the predicted 
velocity structures given the presence of the main fila- 
ment. 

The bottom row of Figure shows the spectra over 
the full density range, i.e. all material contributes to the 
emission. Starting with the case n c = 10 4 cm -3 (bottom 
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Fig. 6. — Gas evolution along ballistic trajectories, for the three 
central filament densities n c = 10 4 , 3 x 10 4 , 10 5 cm -3 as indicated. 
(a) Density compression ratios along all trajectories of Fig. [5] (b) 
Parcel velocities along all trajectories. 

left), we first note that the lines-of-sight #6 and #7 are 
just for demonstration: the spectra are perfectly sym- 
metrical around vi os = by construction. Line-of-sight 
#5 starts out along the direction of maximum compres- 
sion within the R~ 2 regime (see discussion in §3.2.2jh 
intersects the filament, and then transverses a region of 
moderate acceleration (within the R~ x regime). This is 
mirrored in the spectrum: the positive velocities peak 
around vi os = 1.5 km s _1 , while the negative velocities 
show a peak at 0.9 km s _1 . Line-of-sight #4 sees a com- 
pression at lower negative velocities (corresponding to 
the far side), while at positive velocities, it crosses a re- 
gion of stronger compression and stronger acceleration. 

Lines-of-sight #1 and #2 offer comparative diagnos- 
tics: they run through the domain in parallel, with #1 
grazing the filament, and #2 intersecting. The spectrum 
along #1 has a velocity range narrower by a factor of 2 
over spectrum #2, and overall lower emission. Line-of- 
sight #2 is far enough in the filament that it picks up 
the emission from the filament itself, and that it mainly 
tracks the modest density compression gas that falls onto 
the filament. Yet, given its orientation, the radial veloc- 
ity forms a small angle with the infall velocity. For #1, 
because of the curvature of the trajectories, the angle be- 
tween radial and infall velocity is larger, thus, the veloc- 
ity spread is smaller. For a more comprehensive picture, 
position- velocity plots will be useful (S I3.2.4)) . 

With increasing central filament density (center and 
right column of Fig. [7]), the total emission is being dom- 
inated by the filament, and dynamical signatures are in- 
discernible for n c = 10 5 cm -3 . Yet, the assumption of 
one tracer covering the full density range is unrealistic. 
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Fig. 7. — Spectra for two density ranges (labeled 12 CO and C ls O; see text), and for the full density range ("total"), calculated along 



lines-of-sight shown in Fig. [8] Columns correspond to the central filament densities of n c 
intensity units are arbitrary and are given in percent relative to the maximum intensity. 



10 4 ,3 x 10 4 ,10 5 cm -3 as indicated. The 




Fig. 8. — Locations of lines-of-sight used for the spectra shown 
in Figure 

The upper rows in Figure show the spectra for mate- 
rial emitting within specific density ranges. For our pa- 
rameter regime, the gas traced by the low-density tracer 
shows the strongest kinematic differences (top row) . The 
central filament has been nearly completely suppressed, 
leaving the infall signatures visible. With increasing fila- 
ment densities, the infall peaks (which are located closest 
to the filament) appear at increasingly higher velocities. 
Higher-density tracers (second row) start to pick up the 
central filament again. 

3.2.4. Centroid Velocity Signatures 

Figure [9] provides a qualitative picture of the gas kine- 
matics around a filament in free-fall accretion. Each 
panel corresponds to a position-velocity (PV) plot for 
lines-of-sight parallel to those indicated in Figure [51 The 
intensity units are arbitrary, and they result from weights 
proportional to the local density enhancements. Lines- 
of-sight with slopes flatter than 1 in Figure [8] have been 



sampled along the ?/-axis, and those steeper along the x- 
axis, resulting in different "position" -ranges in Figure [9] 
Also, position ranges are limited by the presence of a 
tracer, thus, spatial scales vary. Note that the plane of 
sky in these PV plots is (in Fig. [8]) perpendicular to the 
line-of-sight and the filament. Small-scale noise in the 
maps is due to unresolved velocity gradients and can be 
ignored. Interpolating between cells with velocity dif- 
ferences larger than the thermal width can remove these 
artifacts. Likewise, detailed structures within the PV en- 
velope should not be over-interpreted because of the very 
basic assumptions regarding the density-dependence of 
the tracers. 

Lines #1 and #3 (Fig. [9]) show the density enhance- 
ments in the low-density tracers of the near and far in- 
falling material, at relatively high velocities, i.e. close 
to the filament proper. At position= —4 pc (line #1), 
material close to the left end of the filament is traced, 
while at position= pc, a contribution from near mate- 
rial due to cylindrical infall (minor compression) and far 
material {vi os < 0) towards the right end of the filament 
(thus compressed) is being picked up. Line #5 is an ex- 
treme case, showing the strong density enhancement of 
the near material, while the far side falls in at lower den- 
sities and velocities (especially in projection). The PV 
maps for lines-of-sight #6 and #7 show the symmetrical 
pattern to be expected for the infall velocities projected 
on the line-of-sight. 

Higher-density tracers tend to map gas closer to the 
filament, and thus at higher velocities. For increasing 
central filament densities (not shown), the filament itself 
becomes more prominent. 

4. DISCUSSION 
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4.1. Gas Accretion onto Filaments 

^demonstrates that a complete picture of filamentary 
molecular cloud evolution needs to consider the environ- 
ment of these clouds. Given their masses, they generate 
gravitational potential wells deep enough that it seems 
unrealistic to neglect accretion. For typical cloud param- 
eters, the cloud mass can double within a Myr, given the 
availability of a mass reservoir. 

^erences to Fischera & Martin (2012) 



4.1.1. 

In thei r analytical study of iso thermal hydrostatic 
cylinders, iFischera fc Martini (J2012h develop a model of 
filamentary molecular clouds pressurized by the ambient 
medium. They focus on the structure of the filament 
itself, depending on the external pressure. Instead of im- 
posing a density profile by varying the exponent p, they 
argue that reducing the external pressure will lead to a 
steepening of the profile (see their Fig. 3), from p < 4 to 
p = 4 for zero (vacuum) pressure. Though the profiles in- 
deed flatten with increasing external pressure, the value 
of p - and thus any power-law fit - is actually unaffected 
(Heitsch, in preparation). 

The very assumption that the filaments are embedded 
in a parental molecular cloud mandates not only that 
the density outside the filament is non-zero, but also, 
that there is a (possibly s ubstan tial) external pressure 
component. iKainulainen et all (|2011f ) discuss evidence 
for dense structures in molecular clouds to be pressure- 
confined by the ambient cloud mass - by themselves, 
the cores would be unbound. Yet, here, we are con- 
cerned with the accretion behavior of the filament, and 
less with the filament's radial structure, and thus leav- 
ing the discussion of whether equilibrium considerations 
are appropriate for another time. Instead of deriving 
the filament structure from first principles, we assume 
that we can subsume any physical processes leading to 
a flattening of the profile into setting a profile index, 
consistent with observed values. Again, since the profile 
is only used to determine the line mass m, this "mor- 
phological" approach seems acceptable - our results re- 
garding filament growth through accretion do not change 
qualitatively with different values of p. Flatter than 
isothermal (p = 4) profiles have be en ascribed to mag- 
netic fields (IFiege fc Pudri tz 2000 a]), accretio n, and non- 
isothermality (Naka mura fc Umemural 1999, Heitsch, in 
preparation). The advantage of imposing p is that the 
model should be applicable to a variety of physical pro- 
cesses resulting in different values of p. Thus, the model 
is to some extent agnostic on contentious issues like the 
admissability of hydrostatic equilibrium considerations 
in a turbulent environment. Obviously, imposing p is also 
a disadvantage - the physical properties of the filament 
together with its environment are not fully explored. 

4.1.2. Accretion and "Nessie" 

I Jackson et al.l (j2Q10h estimate typical core separations 
of 4.5pc for the "Nessie" -filament [(/,&) = (338.4,-0.4)] 
of a total length (at a distance of 3.1kpc) of 80pc. As- 
suming an isothermal cylinder, they conclude that the 
core separation is by a factor of « 5 larger than the 
estimated central density of 10 4 cm -3 would warrant. 
The problem is directly obvious from Figure [Tfc and [T]i. 
At a central density of n = 10 4 cm~ 3 , the cor e separa- 
tion should be ~ 1 pc, but I Jackson et all ((2010) observe 



4.5pc, corresponding to a much lower central density of 
(i n our model) 200 cm ~ 3 . 

I Jackson et al.l (J2010h suggest addition al support mech- 
anisms, such as magnetic fields (jFiege fc Pudritzl 
2000a. b), and/or turbulence. As can be seen from 
Figure [3fc and [3]i, only the strong magnetic field case 
(/? = 0.3) has a chance of reaching an equilibrium solu- 
tion (with zero infall). Yet, the densities and the frag- 
mentation length still fall short by a factor of ~ 2. 

In the context of the discussion in §2.3| "Nessie" 
would still be accr eting gas, possibly dr i ving turbu- 
lence via accretion (jKlessen fc Hennebelld l2010f ) . Yet 
while accretion- driven turbulence can slow down accre- 
tion and thus can delay the onset of gravitational frag- 
mentation and collapse, it cannot prevent it (Fig. Q]). 
Even when starting with a turbulent velocity dispersion 
(at rather generous ef ficiencies of 10%, see discussion by 
iKlessen fc Henneb elle 201QJ), the accretion velocity is al- 
ways larger and grows faster than the turbulent veloc- 
ity dispersion. This short-coming of a turbulent support 
scenario is independent of the fact that the turbulent ve- 
locity field itself - in order to give rise to a supporting 
pressure - would have to be uniquely configure d, which i s 
extremely unlikely to occur (jBurkert fc Hartmannl l2004). 
We conclude that while accretion will drive turbulent mo- 
tions in a filament, and while it possibly will lead to frag- 
mentation of the filament, it cannot lead to turbulence 
that in turn can support the filament. 

" Nessie" could have fragmented at an earlier stage dur- 
ing its formation, with the sub-filaments still accreting 
gas and thus growing. To reach the characteristic frag- 
mentation scale of 4.5 pc, a central density of 200 cm -3 
would be needed, i.e. the fragmentation would have had 
to occur very early during the formation in deed. 

Finally, Herschel studies of filaments (e.g. lAndre et al.l 
120101 : lArzoumanian et al.l I20TTI : iPalmeirim et al.l I2013D 
suggest that with higher resolution and sensitivity avail- 
able, more cores and fragmentation centers seem to ap- 
pear, possibly indicating that the 4.5 pc quoted by 
I Jackson et al.l (120101 ) may be an upper limit for the frag- 
mentation scale. 

4.1.3. Filament Evolution and FWHM(N C )- Correlations 

In their Herschel study of dust filaments i n Aquila, 
Polaris and IC5146. rArzoumanian et all (J2011I ) conclude 
from a radial profile analysis that their identified fil- 
aments may share a similar characteristic width of ~ 
0.1 pc. They also point out that the filament width 
does not depend on the central column density N c , even 
for the gravitationally unstable filaments in their sample. 
They argue that a tu rbulent filament form ation mecha- 
nism as discussed by Pa doan et al.l (J2001) may explain 
the similar widths for gravitationally stable filaments, 
and that for gravitationally unstable filaments it could be 
a consequence of continuing accretion, assuming a virial- 
ized filament. 

To test the latter statement, 2800 filament accretion 
models (see $2]) were run, varying four parameters: the 
profile exponent 1.25 < p < 3 (see eq. [TT]h the driving 
efficiency 0.5% < e < 5% (see eq. [T6|) . the background 
density 10 2 < no < 3 x 10 3 cm -3 , and the ratio between 
the initial central density and the background density, 
1.1 < n c (0)/no < 5.0. The choice of ranges is to some 
extent dictated by getting stable solutions. The physics 
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were restricted to accretion-driven turbulence, since the 
isothermal and the magnetized cases do not leave any 
option to break the FWHM(7V C ) correlation. As a cor- 
rollary, it is unlikely that the infall is magnetically domi- 
nated, at least not at the s implistic level considered here 
(also, see Palmei rim et al.ll2Q13L for a discussion of mag- 
netically dominated accretion in Taurus). 

Figure [10] summarizes the test results. In panels (a) 
and (b), each data point represents the characteristic 
values of the FWHM and central column density for one 
model (thus, there are 2800 data points). In Figure [TOk. 
the values are picked at the end of the model, i.e. for 
the final (maximum) column density. Thus, nearly all 
filaments are gravitationally unstable (towards large N c 
with respect to the Jeans length, indicated by the dashed 
line), since the integration proceeded beyond m/m cr = 1. 
For Figure [TUb , the data points are randomly sampled, 
weighed by the time spent in a given column density 
bin. The positions in the (FWHM, 7V c )-plane are deter- 
mined by the median of the column density range within 
which the filament spends 90% of its time. Since the col- 
umn density tends to grow more slowly at early stages, 
the early stages of the filament's growth are sampled. 
Thus, there is a good fraction of gravitationally stable 
filaments. Figure [TOb maps the probability to find a fil- 
ament at a given position in (FWHM,iV c ) space during 
its evolution. The evolutionary tracks of single filaments 
can be guessed from the diagonal streaks towards the in- 
stability line, and parallel to it for large column densities. 

The results are analyzed with the expressions f or the 
(column) density profile by lArzoumanian et all ([2011, 
their eq. 1). Since the exponent p is one of the random 
parameters, we determine the constant A p for a range of 
exponents and fit it with a cubic. The FWHM is deter- 
mined by 



FWHM 



R (V^" 1 ) - l) 



1/2 



(25) 



If we choose to "observe" the model filament at its fi- 
nal stage (Fig. [TOk). the upper envelope of the column 
density distribution is flat. There is a tendency towards 
lower FWHM with higher column densities, mirroring 
the expected effect of filament contraction. Most of the 
filaments' FWsHM reside within the observed range (dot- 
ted lines). Nearly all data points are located in the un- 
stable regime - only a small fraction of filaments never 
makes it to gravitational collapse. 

The time-sampled distribution (Fig. HUb) extends into 
the stable regime, and it suggests consistency with the 
observed decorrelation between FWHM and N c . Clearly, 
the initial (gravitationally stable) conditions of the fil- 
aments are now contributing. The FWHM decreasing 
with smaller column densities is to a large extent an ar- 
tifact of the initial conditions: we choose a random com- 
bination of parameters and assume that the filaments 
have evolved profiles. In other words, these models can- 
not address the question of filament formation. As in 
Figure [TOk . there is a tendency towards smaller FWHM 
at high column densities. 

Finally, Figure [T0h gives a summary view of the fila- 
ment distribution and evolution. As for the time-sampled 
case of panel (b), the filament trajectories in the sta- 
ble regime (left of the dashed line) are possibly a re- 
sult of our initial conditions. Yet, over a large range of 



observed column densities, the FWHM stays flat. We 
overplotted the FWHM( ^) values for selec ted filaments 
in IC 5146, drawn from lArzoumanian et al.l (J2011L their 
Table 1). Colors indicate whether the filament contains 
YSOs (red), pre-stellar cores (blue), cores (green), or 
nothing (white). 
We note several points: 

(a) — Overall, the models can reproduce a flat 
FWHM(iV c ) distribution. The flatness is a direct conse- 
quence of accretion-driven turbulence. For low efficien- 
cies e (blue colors in Figs. ITOk.b) a down-turn at high 
column densities is more readily observed. The isother- 
mal case not including turbulence (not shown) results in 
a strong correlati on between FWHM and N c similar to 
that observed by iFischera fc Martini (J2012h . The same 
effect can be seen for the magnetic cases (see Fig. [2j [3]): 
either, the scaling of the field with density is too weak to 
affect the results (and thus yielding a strong correlation 
between FWHM and N c ), or the filament widths stay 
flat with iV c , but are an order of magnitude larger than 
observed. 

(b) — At high column densities, the FWsHM decrease, 
as one would expect for gravitationally contracting fil- 
ame nts. Yet this is in consistent with Figure 7 of 
lArzoumanian et "all (J2011I ). At such high column den- 
sities, the filaments are expected to fragment and thus 
form stars, resulting in structures that may not be rec- 
ognized as filaments. Thus, if there were a down-turn, it 
may be hard to observe. The following issue is also puz- 
zling and may warrant further discussion: Ignoring the 
filaments with the two highest column densities (log N c = 
22.1,22.2) results in a FWHM(iV c ) distribution that is 
clearly constrained by the Jeans length. In fact, the fila- 
ments containing pre-stellar cores (blue) and YSOs (red) 
nearly all follow the Jeans length, while the remaining 
cores reside well to the left of the Jeans length. As for 
a justification of ignoring t hose two data points , com - 
paring Figure 3b and 4a of lArzoumanian et al.l (|2011h 
suggests that the filaments in question (#6 and #12) 
are in the region with hig hest confusion and th us present 
possibilities for overlaps. lJuvela et al.l (|2012h point out 
that overlaps and three-dimensional effects can affect the 
filament parameters derived from two-dimensional pro- 
jections. As those authors show, what appears to be a 
continous filament in PPV space, may well be a set of 
(isolated) str uctures in three dim ensions (see also Figs. 
2, 3 and 4 of lHeitsch efani2009al ). 

Obviously, in Figure [TUb . we have restricted ourselves 
to the readily availab le data set of IC 51 46, leaving out 
the data for Aquila (jKonyves et al.l 120101 ). on which Ar- 
zoumanians et al. conclusions largely rest. 

(c)— The (FWHM, 7V c )-plane (Fig. [TOb) is most densely 
populated at FWsHM somewhat larger than observed. 
The "center of mass" of the model filaments resides 
more at 0.3 pc, ra t her th an at 0.1 pc, as claimed by 
lArzoumanian et all ([2011). One possible reason could 
be the choice of the driving efficiency e: it is mostly the 
high efficiencies that contribute to the large FWsHM. 
Tests like these could be used to constrain observation- 
ally the turbulent efficiencies, but this would require a 
more detailed (and more realistic) approach to modeling 
the effect of accretion on the filament. The difference 
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between modeled and observed FWsHM could also stem 
from how the FWsHM are measured. Finally, they could 
arise from observational selection effects. Though it is 
tempting to speculate whether observations tend to pick 
out filaments at their late stage of evolution (Fig. ITOk ) 
because of the highest available contrast, or at earlier 
times because of more time being spent at lower column 
densities, the current study cannot give a definite answer. 
Yet, it seems likely that the observed filament samples 
are not unbiased. 

In summary, the filament models including accretion- 
driven turbulence can reproduce the observed decorre- 
lation between filament FWHM and column density. 
For higher column densities, a correlation eventually 
emerges, yet this stage would also be affected by frag- 
mentation and ensuing disruption. For low column den- 
sities, the formation mechanism may play a role ( §4.1.5)) . 

4.1.4. Validity of the Free-Fail Assumption 

The discussion hinges on the assumption that the 
gas is in free-fall towards the filament. This may - 
at first sight - seem an extreme assumption. Yet, for 
the sake of offering an alternative picture to the clas- 
sical equilibrium considerations, the free-fall assump- 
tion is meant as a counter-point. While there is ten- 
tative observational evid ence for accretion onto molec 
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2009), see discussion 
the situation is some- 



ula r clouds (e.g. 
in iKlessen fc Hennebelld 1201 

what less clear in filaments, although recent work 
seems to sup port the cla im tha t there are glob al in- 
fall motions (iMotte et al.l 120051: iSchneider et~aTl 120101: 
iHennemann et all 120121 : iPalmeirim et al.l 12013. see also 
Friesen et al., in preparation). Though somewhat smaller 
than predicted here, the observed kinematic signatures 
are on the order of the virial velocities, which would make 
them mass-dependent. As long as the infall velocities de- 
pend on mass, there is a basis to assume at least some 
free-fall component. 

4.1.5. Initial Conditions 

The choice of the initial central filament density of 
n c (0) = 200 cm -3 is motivated by the goal to high- 
light the evolution of the filament over a large dynamical 
range. As is clear from e.g. Figure [TJd, at such low densi- 
ties, the accretion timescales are comparable to turbulent 
or dynamical timescale, questioning whether considering 
the filament as unaffected by its environment is realis- 
tic. Yet, choosing a higher n c (0) jus t results in a la rger 
= m/m cr as starting point. As iFischera fc Martini 
2012) point out, / can be used as a proxy for the time 
in the filament's evolution. Forming filaments of a few 
hundreds to thousands cm -3 poses no problem in a tur- 
bulent ISM, especially given ( transient) converging flows 
and strong radiative losses (jAudit fc Hennebelle 2005; 



Heitsch et al. 2006; Vazquez- Semadeni et al. 2006, 2007; 



Heitsch et al.l 2008; C arroll-Nellenback et a l. 2013), even 



in the presence of magn etic fields (jlnoue fc Inutsukal 
2008; He itsch et al.ll20 09b), or, in a more realistic set ting, 
by wind-driven superbubbles (Ntorm ousi et al.1l2011f ). 

4.2. Effect of tidal forces on filament accretion 

Tidal forces around finite cylinders can lead to stretch- 
ing and compression of infalling fluid parcels (Q. Inter- 
preting the fluid parcels as "clumpy" molecular gas in the 



filament's vicinity, this effect could explain the "gravita- 
tional streamers" or "f ans" around star- forming filaments 
(see[Myers 2( KMl2QllL for a selection of objects, and also 
for an alternative explanation involving instabilites). 

4.2.1. Observational Tests 

The discussion in ^31 specifically Figure [9j could be 
used in two ways to interpret observations of a filament 
of known mass. First, if one assumed that the gas around 
an observed filament is in approximate free-fall, then the 
variation in the spectra with line-of-sight positions (i.e. 
the PV plots) could be used to determine the orientation 
of the filament with respect to the line-of-sight. Con- 
versely, if the orientation were known, one could esti- 
mate to what extent the gas is actually in free-fall around 
the filament. At minimum, infall signatures could be 
searched for, independent of whether this infall is in free- 
fall or not. Uncertainties are introduced by an assumed 
starting position of the fluid parcels - it is unlikely that 
they will be at rest. 

Figure [9] offer an observable test of the hypothesis that 
the "fans" are caused by tidal forces. Yet, the analysis 
assumes that the line of sight and the filament are located 
in the same plane, i.e. that a projection of the line-of- 
sight has a parallel component to the filament. Thus, for 
a comparison of Figure [9] with observations, the observa- 
tional spectra should be taken along the axis of the cen- 
tral filament (see Friesen et al., in preparation). A more 
detailed study would include full three-dimensional ef- 
fects, with line-of-sight components perpendicular to the 
filament axis. This would also catch the signature of the 
cylindrical infall, if present. This complication has been 
omitted for the sake of clarity, but it should be addressed 
i n the future t o ena ble a more robust observational test. 

iKirk et all (|2Q13h find signatures of flows along fila- 
ments. Such dynamics internal to the central filament 
have been neglected here, simplifying the interpretation 
of the PV plots. Since the infalling material presumably 
shocks at high densities close to the filament, the internal 
motions are likely to be smaller than the infall motions. 
Also not considered was the possibility that the tracers 
could very well be optically thick, although in view of the 
large velocity gradients, this might be less of an issue. 

Given the free-fall velocities of a few km s _1 at the po- 
sition of the filament, the gas could possibly shock. This 
assumes a well-defined filament, and a more or less co- 
herent infall. From the observed (column) density struc- 
tures, it is not clear how realistic such an assumption, 
is, but in any case , it could be test ed with appropriate 
shock models (e.g. iPon et al.ll2012af ). 

4.2.2. Validity of Ballistics Assumption 

The fluid parcels used to determine the amount 
of stretching and compression have been treated as 
pressure-less, i.e. they can suffer infinite compression. 
Removing this simplification will lead to two effects. 
First, the density compression along a trajectory would 
be expected to be more gradual, instead of the strong 
compression observed in Figure [5] Second, as the fluid 
parcels travel towards the central axis of the filament, 
density and pressure will increase, thus decelerating the 
gas eventually. This would lead to a compression along 
a direction perpendicular to the filament axis, and thus 



Gravitational Infall onto Filaments 



13 








j 


1 ' ■ ' i \' 


S \ 


(a) max(N c ) 










\ 


. 


o 

o 


-1 

-2 






<%> 4 










v ' ou <?*f%$>- 












\ °- 










Cturb [ % ] 


- 




0.5 


2.7 


5.0 - 




-3 


i 




, , I , , , , I , , , 


, i , , , , i , , , , r 



cm 

° -2 



-3t 



(b) 90% time 




, [%] 







-3t 



log probability 



20 



21 

log 



22 



23 



20 



21 
log 



22 



23 



20 



21 22 

log N c [cm" 2 ] 



23 



Fig. 10. — FWHM against column density for all model filaments, (a) Quantities determined at final stage (i.e. maximum column 
density), (b) time-sampled quantities, and (c) probability map. The plots can be directly compared to Figure 7 of Arzoumanian et al. 
(|201lD . The range of observe d FW sHM is given by the dotted lines, and the dashed line indicate the Jeans length, Xj = c* / (G/jlttih N c ) . 
The driving efficiency e (see £12.31 and eq. I16J) is indicated by colors. In panel (c), observed filament values are plotted as diamonds, with 
colors indicating whether the fila ment contain s YSO s (red), pre-stellar cores (blue), cores (green), or nothing (white). Error bars indicate 
the uncertainties given by Arzoumanian et al. (2011). 



to a flattening of the parcels, thus enhancing the appear- 
ance of fan-like structures. 



5. CONCLUSIONS 

In an attempt to elucidate the evolution of a filamen- 
tary molecular cloud including its environment, the char- 
acteristic fragmentation and accretion timescales have 
been evaluated under the assumption of free-fall accre- 
tion. Though certainly an extreme view, it is supported 
by observational evidence of virial-like infall of gas in 
massive star- formation regions. Moreover, the filamen- 
tary structure of molecular clouds suggests pressure-less 
collapse to be dominant at least during the formation of 
a cloud, a notion that in view of the thermal properties 
of the gas is not completely outrageoufl 

The analysis presented in $2] and Figures [T]|3] demon- 
strates that the effects of accretion will play a role for any 
molecular cloud with physically reasonable parameters. 
Magnetic fields and/or turbulence can slow down the ac- 
cretion and thus delay the onset of gravitational fragmen- 
tation, but they cannot prevent it. Including the effects 
of accretion-driven turbulence slows down the build-up 
of the filament by approximately a factor of 2. 



Observational evidence for a decorrelation be- 
tween filament FWHM and peak column density 
(jArzoumanian et al.ll201ll ) can be explained by contin- 
uing accretion - as those authors had speculated. 

Tidal forces around finite filaments can lead to stretch- 
ing and compression of clumpy molecular gas being ac- 
creted onto the filament. The resulting gravitational 
streamers are reminiscent of "fans", offering a possible 
explanation of such structures around star formation re- 
gions. Position- velocity plots for such an accretion pro- 
cess offer the opportunity for a comparison with obser- 
vations of massive filaments, either determining the ori- 
entation of the filament with respect to the line-of-sight 
(if free-fall is assumed), or - if the orientation is known 
- testing the hypothesis of free-fall accretion. 

The referee's report was extremely helpful - thank you 
very much indeed to the anonymous referee. I gratefully 
acknowledge support by NSF grant AST 0807305, and 
by NASA NHSC 1008. Insightful comments by and dis- 
cussions with H. Kirk, D. Johnstone and L. Hartmann 
are very much appreciated. I am indebted to S.E. Ragan 
for explaining observational subtleties and for a thorough 
reading of the manuscript. 
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